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Abstract 

p-tubulins are structural components of microtubules and the targets of benzimidazole fungicides used to control 
many diseases of agricultural importance. Intron polymorphisms in the intron-rich genes of these proteins have been 
used in phylogeographic investigations of phytopathogenic fungi. In this work, we sequenced 2764 nucleotides of the 
|3-tubulin gene (Pp tubB) in samples of Phakopsora pachyrhizi collected from seven soybean fields in Brazil. Pp tubB 
contained an open reading frame of 1341 nucleotides, including nine exons and eight introns. Exon length varied 
from 1 4 to 880 nucleotides, whereas intron length varied from 76 to 1 02 nucleotides. The presence of only four poly- 
morphic sites limited the usefulness of Pp tubB for phylogeographic studies in P. pachyrhizi. The gene structures of 
Pp tubB and orthologous p-tubulin genes of Melampsora //n/and Uromyces viciae-fabaewere highly conserved. The 
amino acid substitutions in p-tubulin proteins associated with the onset of benzimidazole resistance in model organ- 
isms, especially at His 6 , Glu 198 and Phe 200 , were absent from the predicted sequence of the P. pachyrhizi p-tubulin pro- 
tein. 
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Phakopsora pachyrhizi Sydow & Sydow, the causal 
agent of Asian soybean rust, is an obligate, biotrophic plant 
pathogen initially identified in Japan (Henning, 1903) and 
nearly a 100 years later in Brazil (Yorinori and Paiva, 
2002), where it accounts for 30%-75% of soybean losses in 
the field. Asian soybean rust has been managed with the use 
of strobilurin, triazole and benzimidazole fungicides (Soa- 
res et al, 2004). However, alternatives are necessary be- 
cause the use of fungicides increases production cost, is en- 
vironment unfriendly, and may trigger pathogen resistance. 
The challenges are enormous as no geographic region or 
state is free of Asian soybean rust and no resistant variety 
has yet been released. In previous work, we reported evi- 
dence that supported an African origin for P. pachyrhizi 
found in Brazil and suggested that multiple, independent, 
long-distance dispersal events were a plausible mechanism 
of introduction (Freire et al, 2008). For a species that was 
introduced recently and that appears to be maintained asex- 
ually (Yorinori and Paiva, 2002), the internal transcribed 
spacer regions (ITS) contain a surprisingly high level of 
intraspecific sequence variation (Freire et al, 2008). To re- 
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fine our inferences on the geographic patterns of genetic 
variation in Brazilian P. pachyrhizi complementary analy- 
ses with a less polymorphic region are necessary. 

As the major components of microtubules, a-tubulins 
and P-tubulins have important cellular functions and are 
amongst the most highly conserved eukaryotic proteins 
(Wade, 2007). Moreover, P-tubulins are a target for benzi- 
midazole fungicides (Davidse, 1986). The mechanism of 
resistance to benzimidazole compounds has been associ- 
ated with specific amino acid changes in the P-tubulin of 
organisms such as Aspergillus nidulans and Neurospora 
crassa (Jung et al, 1992; Cooley and Caten, 1993; Oakley, 
2004). Despite the high level of conservation at the protein 
level, the P-tubulin genes are intron-rich and are therefore 
potentially useful in fungal phylogeography (O'Donnell et 
al., 2000; Ceresini et al, 2007; Kauserud and Shalchian- 
Tabrizi, 2007; Konrad et al, 2007). 

In this report, we (1) describe and analyze the com- 
plete sequence of the P-tubulin gene of P. pachyrhizi (Pp 
tubB), (2) compare the gene structure of Pp tubB and the 
orthologous p-tubulin genes of the rust pathogens 
Melampsora lini and Uromyces viciae-fabae, (3) discuss 
the amino acid sequence alignments and phylogenetic anal- 
yses of P-tubulins of P. pachyrhizi, other rust pathogens, 
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and model organisms, and (4) identify specific amino acid 
residues of (3-tubulin potentially associated with benzi- 
midazole resistance in P. pachyrhizi. 

By using the gene sequence of the M. lini P-tubulin 
gene (GenBank: AF3 17682) (Ayliffe et al, 2001) as the 
initial query, we ran a BLAST search (Altschul et al. , 1 990) 
against public databases at the National Center for Biotech- 
nology Information. The search retrieved several highly 
similar sequences (E-values = 0.0) of filamentous fungi, es- 
pecially basidiomycetes. Among the sequences retrieved 
from GenBank was a 40,837 nucleotide genomic clone of 
P. pachyrhizi. This clone (JGIAFNA-573F20a; GenBank 
AC 170158) contained a segment with 84% identity to the 
query sequence. Further refinement of the sequence align- 
ments showed that this segment contained sequence fea- 
tures which were highly conserved among P-tubulin genes 
of basidiomycetes (data not shown). This segment was fed 
into the program Primer3 (Rozen and Skaletsky, 2000), and 
specific primers were designed for PCR amplification of 
the entire coding sequence, intervenient sequences, and ad- 
jacent flanking regions at the 5' and 3' termini of the Pp 
tubB. With the assistance of the program PrimerSelect 
(DNASTAR, Madison, WI, EUA), ten primer sequences 
were selected (1: 5 '-GACACGGTAAGGGCTTGAGT-3 ' ; 
2: 5'-CAAGGTGCTTCCCACATACC-3'; 3: 5'-CCTCC 
AAAGTGTCAGTCAAAC-3 ' ; 4: 5 '-CGGGGTACATAC 
TTGTTGGC-3'; 5: 5 '-CTCGATCAGGTACAAGGGAA 
C-3'; 6: 5 ' -TC AAAC ATCTGGGAGGTCAG-3 ' ; 7: 5'- 
CCCTACAATGCGGTTCTCTC-3 ' ; 8: 5'-GTTGGACTC 
AGCCTCTGTGA-3 ' ; 9: 5 '-GTGGCAGCTTATTTCAG 
GGG-3'; 10: 5 ' -CC AATTCCCT VTGTTACTGA-3 ' ) . 

Phakopsora-infected soybean leaves were collected 
from seven soybean fields (Table 1). Uredinospores were 
harvested and total genomic DNA was extracted according 
to Freire et al. (2008). PCR amplifications of the Pp tubB 
were done with the primers indicated above, in five combi- 
nations (1-2, 3-4, 5-6, 7-8 and 9-10). Electrophoresis in 
agarose gels indicated that the reactions each yielded a sin- 
gle fragment of the expected size. Non-incorporated nu- 
cleotides and primers were removed using EXOSAP (USB 
Corporation, Cleveland, OH, USA)). 

The primers used in the original PCR amplification 
were also used for sequencing. All of the sequences were 
imported into SEQUENCHER version 4.8 (Gene Codes 
Corp., Ann Arbor, MI, USA) for editing and DNA se- 
quence alignment. A consensus sequence for the entire Pp 
tubB, including flaking regions, was generated for each 
sampling site. The relative positions of the initiation codon, 
exons, introns and termination codon in Pp tubB were in- 
ferred from alignments with published orthologous gene 
copies of other basidiomycetous fungi that were used as 
guides. DNA sequences were translated into protein se- 
quences using the universal code after removing the 
introns. A dataset of protein sequences was constructed us- 
ing the program T-COFFEE (Notredame et al., 2000) with 



Table 1 - Geographical location of the seven soybean fields from which 
Phakopsora-infected soybean leaves were collected. 



Soybean field 


State 


Latitude 


Longitude 


Chapadao do Sul 


Mato Grosso do Sul 


-18°47'39" 


52°37'22" 


Guarda Mor 


Minas Gerais 


-17°46'15" 


47°05'54" 


Palmas 


Tocantins 


-10°12'46" 


48°21'37" 


Piracicaba 


Sao Paulo 


-22°43'31" 


47°38'57" 


Rio Verde 


Goias 


-17°47'53" 


50°55'41" 


Seberi 


Rio Grande do Sul 


-27°28'41" 


53°24'09" 


Vicosa 


Minas Gerais 


-20°45'14" 


42°52'55" 



default values. The phylogenetic relationships of Pp tubB 
protein with P-tubulins of other filamentous fungi, 
Saccharomyces cerevisiae and Glycine max were inferred 
using the neighbor-joining method as implemented in 
MEGA version 4.0 (Tamura et al, 2007). Bootstrap values 
were obtained with 1000 replicates. 

Our strategy for sequencing the Pp tubB yielded five 
overlapping DNA fragments that resulted in a total 
aligned sequence of 2764 nucleotides after assembly. The 
sequence was deposited in GenBank with accession num- 
ber GU354165. An initiation codon (ATG) and a termina- 
tion codon (TAG) delimited an open reading frame of 
1341 nucleotides flanked by segments of 401 and 234 nu- 
cleotides 5' and 3' to the coding region, respectively. The 
promoter region contained the signals characteristically 
involved in the initiation of transcription, such as three pu- 
tative TATA boxes located at - 1 76, -320 and -3 86 and five 
putative CAAT boxes located at -92, -229, -237, -269 and 
-350. A pyrimidine-rich element of 30 nucleotides was 
found in the promoter region between -205 and -175. Such 
a pyrimidine-rich region is commonly associated with po- 
tential transcription factor binding sites in many genes of 
filamentous fungi (Cooley and Caten, 1993). The 
AATAAA recognition signal associated with the efficient 
formation of mRNA 3 ' termini (Wahle and Riiegsegger, 
1999) was present 44 bases downstream from the termina- 
tion codon. 

The Pp tubB gene contained nine exons (14, 19, 16, 
55, 49, 880, 172, 121 and 15 nucleotides long, respectively) 
and eight introns (125, 76, 86, 102, 85, 84, 99 and 128 nu- 
cleotides long, respectively). The consensus sequences 
GT(A/g)NGT at the 5' terminus and (C/T)AG at the 3' ter- 
minus flanked the introns. These consensus sequences are 
found at the intron-exon boundaries of many filamentous 
fungi (Begerow et al. , 2004). The seven Pp tubB sequences 
diverged at only four polymorphic sites. We found a 1 bp 
indel in a polyT at 1580, which was either nine or ten bases 
long. The presence of sequences differing in length resulted 
in peak displacements and chaotic electrophoretograms 
when this polyT was reached during the sequencing of 
some samples. We also found three base substitutions, 
identified as ambiguities: R (G or A) upstream from the ini- 
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tiation codon at position -50, M (A or C) in the second 
intron at 231, and Y (C or T) in the third intron at 323. 
These three polymorphic sites were identified as overlap- 
ping double peaks in electrophoretograms of both strands. 
These double peaks appeared in electrophoretograms of 
DNA sampled from distinct fields, which suggested that 
our samples contained a pool of strains and that distinct 
fields shared these strains. Although the Pp tubB gene con- 
tained a large number of introns and intron length varied 
substantially, the phylogeographic usefulness of these fea- 
tures in P. pachyrhizi was hampered by the extremely low 
levels of intron variability within and among the seven soy- 
bean fields sampled. For some organisms, such as 
Coniophora arida and Coniophora olivacea (Boletales), 
the ITS region is less polymorphic than in P-tubulin genes 
(Kauserud and Shalchian-Tabrizi, 2007). The lack of poly- 
morphic sites in the intron-rich Pp tubB gene contrasted 
sharply with the high level of sequence variation we found 
in the ITS1 and ITS2 regions (Freire et al, 2008). 

The gene structures of Pp tubB and orthologous p-tu- 
bulin genes of M. lini and U. viciae-fabae were highly con- 
served in length, distribution, and intron and exon insertion 
sites (data not shown). Currently, M. lini and U. 
viciae-fabae are the only two rust pathogens for which 
complete P-tubulin gene sequences available in public da- 
tabases. Although intron length and position were con- 
served, intron sequences were highly divergent among 
these three rust species. In most attempts, sequence align- 
ments among introns of equivalent position simply were 
not achievable. For those that did align, the maximum iden- 
tity was < 30%. The number of introns in other fungal 
P-tubulin genes varies considerably. For example, there are 
only three introns in the Septoria nodorum p-tubulin gene 
(Cooley and Caten, 1993) whereas the P-tubulin genes of 
many basidiomycetous fungi contain 6-10 introns (Ayliffe 
et al, 2001). 

The Pp tubB gene encoded a full-length protein of 
447 amino acids (Figure 1). A search at the conserved do- 
main database (Marchler-Bauer et al. , 2007) revealed that 
the P. pachyrhizi P-tubulin protein (Pp TUB) had four do- 
mains characteristically found in other P-tubulins: a nucle- 
otide-binding site, taxol-binding site, an a/p domain inter- 
face and a p/a domain interface. Sequence alignment using 
the BLAST algorithm (Altschul et al, 1990) revealed that 
Pp TUB shared high levels of sequence identity with p-tu- 
bulin from other basidiomycetous fungi: U. viciae-fabae 
(98%), M. lini (96%), and Pleurotus sajor-caju (89%), 
other Ascomycota: N. crassa (82%), A. nidulans (82%), 
and S. cerevisiae (73%), higher animal species, such as Bos 
taurus (85%) and higher plant species, such as G. max 
(40%). The amino acid sequence alignment of p-tubulins 
indicated that sequence heterogeneity was greatest in the 
carboxy-termini. This higher level of sequence divergence 



in the carboxy-termini has been associated with the loca- 
tion of this region at the surface of the heterodimers, its role 
as a binding surface for microtubule-associated proteins, 
and its function as a target for post-translational modifica- 
tions (Wade, 2007). 

The mechanism of resistance to benzimidazole com- 
pounds has been associated with specific amino acid substi- 
tutions in the P-tubulin of many organisms. The benzi- 
midazole sensitivity of microsporidia such as Vittaforma 
corneae and Enterocytozoon bieneusi is related to the pres- 
ence of six amino acids: His 6 , Ala 165 , Phei 6 7, Glui 9g , Phe 2 oo 
and Arg 2 4i (Franzen and Salzberger, 2008). Resistance to 
benzimidazole in a S. nodorum mutant arose as a result of a 
base substitution at the first position of codon 6 (C — > T) of 
the P-tubulin gene, which led to the amino acid substitution 
His 6 — > Tyr 6 (Cooley and Caten, 1993). Likewise, benzi- 
midazole resistance in a N. crassa mutant appeared as a re- 
sult of an amino acid substitution at another site 
(Glui9s — » Glyi9 8 ) (Fujimura et al, 1992). Most interest- 
ingly, experiments reported by Jung et al. (1992) with A. 
nidulans showed that amino acid substitutions at some 
sites, such as Ala 165 , conferred benzimidazole resistance 
but lead to abnormal cell growth, whereas amino acid sub- 
stitutions at three other site (His 6 , Glui 98 and Phe 2 oo) con- 
ferred resistance to benzimidazole agents without causing 
any conditional blockage of cell growth. Given that P-tubu- 
lins are amongst the most highly conserved eukaryotic pro- 
teins, Oakley (2004) fitted the amino acid sequences of A. 
nidulans a- and P-tubulins to the three-dimensional struc- 
ture of Bos taurus a- and P-tubulins (Nogales et al, 1998) 
and showed that amino acids associated with benzimi- 
dazole resistance in A. nidulans (His 6 , Tyr 50 , Gln 134 , Alai 65 , 
Gluigg, Phe 2 oo and Met 2 57) were structurally close together; 
this finding suggested that the binding site for benzimi- 
dazole agents may be in the region defined by these muta- 
tions. 

To the best of our knowledge, no studies have exam- 
ined the relationship between amino acid substitutions 
and benzimidazole resistance in P. pachyrhizi. As shown 
here, the amino acid substitutions in P-tubulin associated 
with benzimidazole resistance in closely related organ- 
isms were absent in the predicted Pp TUB sequences from 
the P. pachyrhizi samples collected in seven Brazilian 
soybean fields. Nevertheless, we found a high level of se- 
quence identity between Pp TUB and P-tubulins of model 
organisms such as A. nidulans and N. crassa (Figures 1 
and 2), for which there is a clear relationship between 
amino acid substitutions at given sites and benzimidazole 
resistance (Jung et al., 1992; Cooley and Caten, 1993; 
Oakley, 2004). We hypothesize that the emergence of re- 
sistance in P. pachyrhizi is likely to follow a similar pat- 
tern of amino acid substitutions in Pp TUB. However, we 
cannot exclude the possibility that there were no 
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Figure 1 - T-Coffee alignment of (3-tubulin proteins. The amino acid residues are shown in single-letter code. Dots indicate residues identical to those of 
Phakopsora pachyrhizi (3-tubulin and hyphens indicate gaps. Conserved domains are indicated as follows: (*) beta/alpha domain interface, (-) alpha/beta 
domain interface, (#) alpha/beta and beta/alpha domain interface, (+) nucleotide binding site and (A) taxol binding site. Black inverted triangles indicate 
amino acids thought to be associated with benzimidazole sensitivity in A. nidulans. GenBank accession numbers are: P. pachyrhizi (GU354165, this 
work), M. lini (Q9HFQ3), U. viciae-fabae (Q96TU8), P. sajor-caju (AAD21093), N. crassa (XP 957669), A. nidulans (XP 658786), S. cerevisiae 
(V01296) and S. taurus (1JFF). 



benzimidazole resistant mutants in the soybean fields at 
the time of sampling. Because of the limitations inherent 
to our sampling and sequencing strategies, we can only 
imply that if mutants were indeed present, their frequen- 
cies were extremely low. The appearance of natural mu- 



tants with amino acid substitutions in Pp TUB, especially 
at His 6 , Gluips and Phe 2 oo, should be closely monitored 
since the continual use of benzimidazoles for controlling 
Asian soybean rust may well increase the selective pres- 
sure under field conditions in the near future. 
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Figure 2 - Neighbor-joining tree based on P-tubulin proteins. Relative 
branch lengths are shown to scale. Bootstrap values (1000 replicates) are 
indicated. GenBank accession numbers are: P. pachyrhizi (GU354165, 
this work), M. lini (Q9HFQ3), U, viciae-fabae (Q96TU8), P. sajor-caju 
(AAD21093), N. crassa (XP 957669), A. nidulans (XP 658786), S. 
cerevisiae (V01296) andB. taunts (1JFF). 



Acknowledgments 

This work was supported by FAPEMIG (grants nos. 
CAG-1484/05 and CAG-0673-3. 08/07 to LOO) and 
IICA/PROCISUR (grant to AMRA). MCMF and TEFM 
were supported by fellowships from FAPEMIG. The au- 
thors thank anonymous reviewers for valuable comments 
on an earlier version of the manuscript. 

References 

Altschul SF, Gish W, Miller W, Myers EW and Lipman D J (1 990) 
Basic local alignment search tool. J Mol Biol 215:403-410. 

Ayliffe MA, Dodds PN and Lawrence GJ (2001) Characterization 
of a |3-tubulin gene from Me lampsora lini and comparison of 
fungal P-tubulin genes. Mycol Res 105:818-826. 

Begerow D, John B and Oberwinkler F (2004) Evolutionary rela- 
tionships among P-tubulin gene sequences of basidiomy- 
cetous fungi. Mycol Res 108:1257-1263. 

Ceresini PC, Shew HD, James TY, Vilgalys RJ and Cubeta MA 
(2007) Phylogeography of the Solanaceae-infecting Basi- 
diomycota fungus Rhizoctonia solani AG-3 based on se- 
quence analysis of two nuclear DNA loci. BMC Evol Biol 
7:163-184. 

Cooley RN and Caten CE (1993) Molecular analysis of the Septo- 
ria nodorum P-tubulin gene and characterization of a beno- 
myl-resistance mutation. Mol Gen Genet 237:58-64. 

Davidse LC (1986) Benzimidazole fungicides: Mechanism of ac- 
tion and biological impact. Annu Rev Phytopathol 24:43-65. 

Franzen C and Salzberger B (2008) Analysis of the P-tubulin gene 
from Vittaforma corneae suggests benzimadozole resis- 
tance. Antimicrob Agents Chemother 52:790-793. 

Freire MCM, Oliveira LO, Almeida AMR, Schuster I, Moreira 
MA, Liebenberg MM and Mienie CMS (2008) Evolutionary 
history of Phakopsora pachyrhizi in Brazil based on nucleo- 
tide sequences of the internal transcribed spacer region of 
the nuclear ribosomal DNA. Genet Mol Biol 3 1 :920-93 1 . 

Fujimura M, Oeda K, Inoue H and Kato F (1992) A single 
amino-acid substitution in the beta-tubulin gene of Neuros- 
pora confers both carbendazim resistance and diethofencarb 
sensitivity. Curr Genet 21:399-404. 



Henning P (1903) Einige neue japanische Uredineen IV. Nova 
HedwigBeih 42:107-108. 

Jung MK, Wilder IB and Oakley BR (1992) Amino acid alter- 
ations in the benA (beta-tubulin) gene of Aspergillus nidu- 
lans that confer benomyl resistance. Cell Motil Cytoske- 
leton 22:170-174. 

Kauserud H and Shalchian-Fabrizi K (2007) Multilocus sequenc- 
ing reveals multiple geographically structured lineages of 
Coniophora arida and C. olivacea (Boletales) in North 
America. Mycologia 99:705-713. 

Konrad H, Stauffer C, Kirisits T and Halmschlager E (2007) 
Phylogeographic variation among isolates of the Sirococcus 
conigenus P group. For Pathol 37:22-39. 

Marchler-Bauer A, Anderson JB, Derbyshire MK, DeWeese- 
Scott C, Gonzales NR, Gwadz M, Hao L, He S, Hurwitz DI, 
Jackson JD, et al. (2007) CDD: A conserved domain data- 
base for interactive domain family analysis. Nucleic Acids 
Res 35:237-240. 

Nogales E, Wolf SG and Downing KH (1998) Structure of the 
aP-tubulin dimer by electron crystallography. Nature 
391:199-203. 

Notredame C, Higgins DG and Heringa J (2000) T-Coffee: A 
novel method for fast and accurate multiple sequence align- 
ment. J Mol Biol 302:205-217. 

O'Donnell K, Kistler HC, Tacke BK and Casper HH (2000) Gene 
genealogies reveal global phylogeographic structure and re- 
productive isolation among lineages of Fusarium 
graminearum, the fungus causing wheat scab. Proc Natl 
Acad Sci USA 97:7905-7910. 

Oakley BR (2004) Tubulins in Aspergillus nidulans. Fungal 
Genet Biol 41:420-427. 

Rozen S and Skaletsky HJ (2000) Primer 3 on the WWW for gen- 
eral users and for biologist programmers. In: Krawetz S and 
Misener S (eds) Bioinformatics Methods and Protocols: 
Methods in Molecular Biology. Humana Press, New York, 
pp 365-386. 

Soares RM, Rubin SAL, Wielewicki AP and Ozelame JG (2004) 
Fungicidas no controle da ferrugem asiatica (Phakopsora 
pachyrhizi) e produtividade da soja. Cienc Rural 34:1245- 
1247. 

Tamura K, Dudley J, Nei M and Kumar S (2007) MEGA4: Molec- 
ular Evolutionary Genetics Analysis software v. 4.0. Mol 
Biol Evol 24:1596-1599. 

Wade RH (2007) Microtubules: An overview. Meth Mol Med 
137:1-16. 

Wahle E and Riiegsegger U (1999) 3 '-End processing of pre- 
mRNA in eukaryotes. FEMS Microbiol Rev 23:277-295. 

Yorinori JT and Paiva WM (2002) Ferrugem da soja: Phakopsora 
pachyrhizi Sydow. Embrapa Soja, Londrina, 6 pp. 

Internet Resources 

Primer 3 software: http://www.genome.wi.mit.edu/cgi-bin/ 

primer/primer3 xgi. 
T-COFFEE software: 

http://www.igs.cnrs-mrs.fr/Fcoffee/tcoffee_cgi/index.cgi. 

Associate Editor: Sergio Olavo Pinto da Costa 

License information: This is an open-access article distributed under the terms of the 
Creative Commons Attribution License, which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



